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Abstract 

We present an analytical method of studying "extended" electronic eigenstates of a diamond 
hierarchical lattice which may be taken as the simplest of the hierarchical models recently 
proposed for stretched polymers. We use intuitive arguments and a renormalization group 
method to determine the distribution of amplitudes of the wavef unctions corresponding to 
some of these "extended" eigenstates. An exact analysis of the end-to-end transmission prop- 
erty of arbitrarily large finite lattices reveals an anomalous behaviour. It is seen that while 
for a special value of the energy the lattice, however large, becomes completely transparent 
to an incoming electron, for the other energy eigenvalues the transmission decreases with 
system size. For one such energy eigenvalue we analytically obtain the precise scaling form 
of the transmission coefficient. The same method can easily be adopted for other energies. 

PACS No.s : 71.23 An, 64.60 Ak, 72.80. Le 



1. Introduction 

Hierarchical lattice models have played an important part in understanding the statistical 
mechanics of phase transitions. Such lattices differ, both in topology and in geometry, from 
the Bravais lattices [1] . One major feature of such structures is their scale invariance, which 
enables an exact implementation of real space renormalization group (RSRG) techniques 
[1-4]. Apart from statistical mechanics, studies of the electronic spectrum of such lattices 
have also revealed striking properties, not shared by common crystalline structures [5-7]. 
For example, Domany et. al. [5] solved the Schrodinger equation on a variety of hierarchical 
lattices using an exact recursive scheme. In the thermodynamic limit, the energy levels were 
found to be discrete, very closely spaced and highly degenerate [5] . The electron localization 
problem however, is not easy to understand on hierarchical lattices. The fluctuating local 
environment around each lattice point is likely to localize the electronic wavefunctions in 
most situations [5] , though for some hierarchical or fractal structures like the Seirpinski gas- 
ket and the Vicsek fractal, "extended" electronic states have been reported [8,9]. This aspect 
makes the study of electronic properties of hierarchical lattices interesting. Normally, the 
"extended" character of electronic wavefunction is associated with translational invariance 
of the underlying lattice and, a hierarchical lattice (like the present one) does not have any 
translational periodicity by virtue of its construction (see Fig. 1). However, we are famil- 
iar with several one dimensional examples, viz. the one-dimensional random dimer [10] or 
quasiperiodic lattice models [11] where one comes across situations in which local positional 
correlation between constituent 'atoms' gives rise to a finite [10] or infinite [11] number of 
"extended" electronic states even though these lattices are not periodic. In hierarchical lat- 
tices, such local positional correlation is not always obvious. On the other hand, the topology 
of the lattice in certain earlier cases has been shown to play an important role in sustaining 
"extended" electronic wavefunctions [8,9]. Such states sometimes exhibit anomalous trans- 
port properties, in the sense that though the amplitude of the wavefunction remains non-zero 
even at distant parts of arbitrarily large lattices, the end-to-end transmission of any incident 
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wavepacket sometimes displays a power law decay [9] . 

Above exact results are available for hierarchical models in which there is a finite variety of 
the nearest neighbour environment around a lattice site (strictly speaking, in any hierarchical 
lattice all sites are inequivalent if one looks beyond the nearest neighbours). The situation 
may thus become quite challenging if one tries to explore a case where the the range of the 
coordination number of the lattice points increases with the generation of bigger and bigger 
lattices. One such example is the well known diamond hierarchical lattice [1] in which the 
coordination number of the vertices range from 2; = 2 to z^ax = 2^"-*^ for an N*^ generation 
lattice. The thermodynamics of spin models on a diamond lattice has been studied exactly by 
RSRG methods [1] . But its electronic properties has not been really studied in detail, a part 
of which we intend to address in the present communication. Our interest in the diamond 
lattice is two-fold. Firstly, we wish to investigate analytically, if there exists any "extended" 
electronic eigenstate though there is no translational periodicity in this lattice. Moreover, 
if such states do exist then, to our mind, it would be interesting to study their amplitude 
profiles and also to study the transmission properties of arbitrarily large finite lattices at 
those special values of the electron-energy for which the distribution of the amplitudes is non- 
trivial. Secondly, we note that similar hierarchical structures have recently been proposed 
by Samukhin et. al. [12] as a possible basic structure of stretched polymers. In stretched 
polyacetylene the polymer network is constructed by coupled polymer chains oriented along 
some direction [12]. The {m,n) hierarchical pseudolattice that serves as the model for such 
polymers may be composed, according to Ref. [12], by taking n bonds forming a chain of 
{n + 1) atoms and then joining m such chains in parallel. It is then quite obvious that our 
diamond hierarchical structure is a (2,2) model lattice in this group (see Fig. 1). So, an 
analytical approach to study the transport properties of a diamond structure seems to be an 
interesting step towards the understanding of the properties of the general (m, n) structure. 

Recently, Zhu et. al. [13], numerically calculated the transmission coefficient of the 
(m, n) polyacetylene model as a function of electron energy for different (m, n) values and at 
different stage numbers. Interestingly, they found, among very fragmented patterns, many 
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energy values for which the transmission coefficient turns out to be unity (or very nearly 
so) for upto fifth generation structures. Also, they have reported (numerically obtained) 
scaling behaviour of conductance for the (m, n) structures, though its precise form is not 
known. In view of this wc consider the simplest (2, 2) version, i.e., the diamond lattice and 
make an analytical attempt to see if such a structure ever becomes completely or partially 
transparent to an incoming electron. Most interestingly, we find that it is possible to get a 
special energy eigenvalue for which arbitrarily large sized diamond hierarchical lattices will 
have transmission coefficient equal to unity. We work out a scheme based on an intuitive 
approach to determine the distribution of the amplitudes of the wave function (ipi) for this 
special energy on a small sized lattice and then use an RSRG approach to work out the 
distribution in lattices for higher generations. The central idea is not too difficult to extend 
to the general (m, n) case. Using the same RSRG formalism we evaluate a whole hierarchy 
of energy eigenvalues for which we again have non-trivial distribution of the amplitudes ipi 
over the entire lattice. The end-to-end transmission across lattices of gradually increasing 
size is now found to decay for this second group of eigenvalues. We explicitly calculate the 
scaling behaviour for one such case and obtain a precise form of the power law followed by 
the transmission coefficient. Scaling forms for other energies can be obtained following the 
scheme prescribed by us. In Sec. 2 we describe our model and the method, in Sec. 3 the 
transmission coefficient is discussed and we draw conclusions in Sec. 4. 

2. The Model and the Method 

In Fig. 1, we show the first three stages of construction of the diamond lattice following 
Berker and Ostlund [2] . A detailed discussion on the construction and topology of lattices 
belonging to this class is given in Ref. [1]. The atomic sites sit at every vertex of the lattice. 
To study the electronic properties of the diamond lattice, we work with the tight binding 
hamiltonian with nearest neighbour approximation: 

H^^e,\i){i\ + Y^U,\z){j\ (1) 

i <ij> 
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where is the on-site potential and tij is the nearest neighbour (n.n) hopping integral. 

As we are primarily interested in the topological aspect of the system we set identical 
values to all the site energies, i.e., = e and all the hopping integrals are taken to be the 
same, i.e. tij = t (chosen to be the scale of energy). However, to facilitate the renormalization 
group calculations we shall designate the site energy of a site with coordination number z, 
as €z. Therefore, in an N*'^ generation lattice e^'s range from €2 to €z^^^, where z^ax — S^""-*^. 

Let us, first of all, try to give a simple intuitive picture of how to construct a wavef unction 
that will have a non-trivial disribution of amplitude over the entire lattice, irrespective of 
its size. We work with the difference equation version of the Schrodinger equation, 

{E-ei)i^i^t (2) 

je{n.n of i) 

where, t/jj is the amplitude of the wave function at the j*'* site. We shall be looking for 
energy eigenvalues that solve the above equation consistently all over the lattice and yet 
yield a non-trivial distribution of iJj/s. We begin with = 3, i.e., where we have at least 
one class of sites with z > 2 { here, Zmax is 4). By inspection we can see that if we choose 
E — €2, then this equation satisfies the Schrodinger equation consistently at all sites if we 
demand that the amplitude of wave function vanishes at sites with coordination number 
four. Extending this idea we find that indeed we can satisfy Schrodinger equation locally 
at all sites in any arbitrarily large lattice provided we set the amplitude ipi at a site with 
coordination number z, equal to zero for all z > 2. That is, in this scheme, the electron 
cannot "feel" the presence of the sites with z > 2. This helps in attaining an "extended" 
character of the wave function. It may be mentioned here that, a similar technique was used 
earher in the case of a Vicsek fractal [9]. Here, it is to be noted that the solutions of the 
Schrodinger equation on such a lattice will be highly degenerate and we are exhibiting one 
such case only. In Figs. 2a and 2b, the amplitude distributions ior E — €2 — are given on 
a lattice with N — 3 and N — A respectively. The pattern of distribution of the amplitudes 
for the N — A case has been obtained by joining the extremeties A and B (Fig. 2a) of the 
basic N — 3 plaquette (which now becomes the building block for the next generation) side 
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by side such that the vertex B of one plaquette falls on the vertex A of the next (Fig. 2b). 
Following this strategy we can use the N — A plaquette as the basic unit and determine the 
amphtude-distribution on an iV = 5 lattice by joining the N = A plaquettes at the vertices 
having z = S. where the amplitude of the wavefunction is zero. The method can be extended 
easily to construct the distribution pattern for higher generations. 

The above idea can now be coupled to an RSRG scheme to extract other possible energy 
eigenvalues. Starting from any arbitrarily large finite version of the hierarchical diamond 
lattice, we can renormalize the lattice n times. The recursion relations for the site energies 
and hopping integral are respectively given by : 

2zt^{n-l) 

e.(») = e.,(„-l)+|^-^L_2_ (3) 

and 

, , 2t2(n- 1) 
t(n) = 7- — S '-— (4 

Here, tzin) and t{n) denote respectively the values of the on-site potential and the hopping 
integral, at the n*'* stage of renormalization. 

Now, suppose we start with an N*^ generation lattice. We renormahze this lattice n — N — 3 
times and bring it down to an "effectively" third generation lattice where, z — 2 and 4. 
The site energies and the hopping integral for this renormalized version are calculated using 
Eqns. (3) and (4). We then apply the earlier trick on this renormalized lattice, i.e., we 
make the amplitudes ipi vanish at the z — 4 sites ( on the rescaled version ). This happens 
for E — £2(71). This is a polynomial equation in E. Though, it is difficult to provide a 
complete proof, yet we have performed explicit calculations upto N — 5 and for n = 1 
and 2. In each case it turns out that the solutions of the equation E = €2(11) satisfy the 
Schrodinger equation all over the lattice with a non-trivial distribution of the amplitudes of 
the wavefunction. Therefore, it is tempting to make the conjecture that the real roots of the 
equation E — 62(71) will correspond to the "extended" eigenstates in the general case, if we 
follow the same prescription. When mapped onto the original lattice the wavefunction will 
vanish at sites with some value of z onwards and will remain finite at all other lattice points 



with lower values of z. Let us clarify this idea by discussing two specific situations. First, we 
choose E — e^iX)- Now we should have at least an = 4 lattice as our starting structure, 
where Zmax = 8. We then renormalize it once to cast it into an effective TV = 3 - stage lattice 
with Zmax = 4. Wc find that the solutions of the equation E = £2(1) are E = ±2, with all 
— initially and t — 1. Each of these energy values consistently satisfies the Schrodinger 
equation everywhere on the renormalized lattice provided we fix ■^j = at the z — A sites 
on this renormalized structure. When mapped onto the original N — A - stage hierarchical 
structure, we see that the amplitude of the wavefunction vanishes only at the vertex with 
the highest coordination number, i.e., 8 and it is non-zero at all other points. The rule for 
constructing such a distribution may be formulated as follows. We consider two plaquettes, 
type / and type //, each being an N — 3 diamond lattice. In Fig. 3a, the plaquette / is 
shown. In type // each non-zero amplitude of the wavefunction is of opposite sign to that 
at the corresponding vertex in type /. In both of them, the Schrodinger equation is satisfied 
at each vertex for E = 2. These two plaquettes, / and //, are then joined end-to-end in an 
alternate fashion such that the vertices with ipi — fall on each other (Fig. 3b). In this way 
the N — A structure is built, and by selecting the N — A structure as the basic unit we can 
construct the pattern for an = 5 lattice by joining them suitably at the extreme points 
(with ijji = 0). The scheme can go on for other higher order lattices. It is to be appreciated 

that, for higher generations ipi will remain zero at all vertices with z = 8, 16, 32, The 

even number of the nearest neighbours helps maintaining the value of ipi equal to zero at 
these sites. For the other sites with 2; < 4 at any generation, the Schrodinger equation is 
locally satisfied everywhere, once it is satisfied at a lower generation. 

Second, we consider E = €2(2). Here, we need to start from an TV = 5 structure. We 
map this lattice onto an TV = 3 version with renormalized parameters and force ipi = 0, as 
before, at the sites with 2; = 4 on this renormalized lattice. The effective TV = 3 lattice 
with the distribution of the amphtudes is presented in Fig. 4a. When unfolded to retrieve 
the original N — 5 lattice, we now find that the amplitude is zero only at the vertices with 
z — 16, and three new amplitudes ±a, ±6 and ±c appear ect z — 2 and 4 - sites. One quarter 
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of the full N — 5 lattice is shown in Fig. 4b with the distribution of the values a, b and c. 
The complementary portion with —a, —b and — c is not shown, but can easily be conceived 
of. In order to satisfy the Schrodingcr equation consistently at each vertex, we see that a, b 
and c should take values {E^ — 8)/8E, {E"^ — 8)/8 and E/8 respectively, provided the energy 
£■ is a solution of the equation E^ — 12E'^ + 16 — 0. This is precisely the polynomial equation 
that is obtained by setting E — £2(2). We thus confirm that for every root of this equation 
we are able to construct extended eigenstates consistent with the Schrodinger equation. We 
expect that, reasoning in the same manner as in the above cases, we are likely to uncover 
a whole set of eigenvalues by solving the equation E = e2{n) with n = 1, 2, 3, in an 
arbitrarily large lattice for which the wavefunction will be "extended" in the sense described 
earlier. 

Before ending this section we point out two other possibilities of getting "extended" type 
of states: 

(i) We consider the third generation lattice, where z„iax = 4. It is quite obvious that by 
choosing E — (which happens to be equal to 62 in the bare length scale in our model), 
we can have a consistent solution of the Schrodinger equation on this lattice in which the 
amplitude of the wavefunction is zero on all vertices with z — 2 and alternates between ±1 
on the vertices with z = A. This is, of course, one of the possible configurations. But we 
can carry on the process of construction for other eigenvalues by setting E = 64,(71) for an 
N — n + 3 generation lattice, and by demanding that ipi vanishes identically on all 2; = 2 
vertices on the n step renormalized version of the same lattice. In Fig. 5, we exhibit the 
distribution of amplitudes for E — 64(1) on a fourth generation lattice. Distribution for 
higher generations and other eigenvalues can be obtained by extending the earlier ideas. 
Similar results can be obtained for the general case with E = e2z{n). 

(a) The other possibility refers to a specific initial choice of the site energies. By looking 
at the recursion relations (3) and (4) for the on-site term and hopping integral respectively, 
we find that if we start with a model where 62^ — ez — zt, then for E — €2 + 2t, each on-site 
term and the hopping integral exhibit a fixed point behaviour. It implies that the nearest 
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neighbour hopping integral does not flow to zero under iteration and we have a non- vanishing 
'connection' between nearest neighbouring sites at aU length scales. This is a clear signature 
of the corresponding cigcnstate being extended [14]. 

We will now describe how to investigate the transmission characteristics of a diamond 
hierarchical lattice. We will emphasize on the analytical treatment of the recursion relations 
and will discuss the behaviour of the transmission coefficient T{E) for the two cases with 
E — €2 and E — 62(1) respectively, for which the transmission coefficient displays totally 
opposite characteristics. T{E) for the other energies {E = e2z{n)) can be obtained by 
following the method adopted for these cases. 

3. Analysis of the Transmission Coefficient 

For calculating the transmission coefficient T{E), we attach two semi-inflnite perfectly 
ordered leads to the two 'diametrically' opposite vertices having the maximum coordination 
number in any generation A^. The original lattice is then renormalized n(= N — 2) times, 
so that we are left with a basic N — 2 rhombus with four vertices each having an effective 
site energy 62(71) and the nearest neighbour hopping integtral t{n) [Fig. 6]. This elementary 
rhombus is now folded into a 'dimer' with site energy and nearest neighbour hopping integral 
respectively given by 

2f2(") 

e = €2 n + 



[E-e2{n)] 
and 

2t^(n) 



[E-e2{n)] 

Following the standard procedure [15] it is then easy to show that 

4 sin^A; 

[P21 - P12 + {P22 - Pn)coskY + (Pii + P22Tsm''k 
where the elements of the matrix P are 
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E-~e^' 



- 1 



t 

To 



(6) 
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and 



^ ) = -^21 (7) 



P22 = -j (8) 



Here, k = cos~^[(ii^ — eo)/2to], where, eo and to refer to the on-site potential and the 
hopping integral respectively, of the ordered lead. 

In order to understand the behaviour of T{E) for large systems at any particular energy, 
we must analyze how the matrix elements Pij behave for large number of RG iterations, 
n. As a lattice of any generation should finally be reduced to a basic N = 2 rhombus 
having only e2{n) and t{n) (see Fig. 5), we must analyse the flow patterns of 62 and t under 
successive RSRG iterations. The evolution of these two parameters ultimately controls e and 
i, and hence the matrix elements Pij. Let us discuss it for two specific cases. Throughout 
the analysis we will set all Sz — and t — 1. 

Case (i) E = €2 = 

From direct calculations we find that as E ^ 0, the leading behaviours (in E) of t and 
62 are given by, 

t(n) ~ ^ (9) 

62(71) ~ (10) 

for n = 2,3,4 and 5. We thus assume these forms to be true for any arbitrary value of 
n, viz., n = m. Then proceeding according to the standard method of induction, we can 
prove, using the recursion relations (3) and (4), that Eqns. (9) and (10) indeed hold good 
for n = m + 1 as well. Therefore, we accept the above forms as the leading terms in the 
expressions for t{n) and €2(71) for — > and for any arbitrary value of n with n > 2. It is 
now easy to work out an expression for {E — e)/i, which is given, to the leading order in E, 

by 

E-i E^ 
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A direct subtitution of the above result in the expressions of shows that, Pu — P22 — 
and P12 — —P21 = — 1 as £■ — > 0. The expression for T{E — 0) now becomes : 

T{E = 0) = ^ '^""'^ , = sin^A; = 1 - A (12) 

If we select eo = 0, then T{E — 0) is unity, and any arbitrarily large diamond hierarchical 
lattice becomes completely transparent to an incoming electron with E — 

Case (ii) E = 62(1) = 2 

The central idea of the analysis for case (i) can now easily be extended to study the 
general situation, where E — £2(71). The results, however, turn out to be totally different, 
as we have checked numerically by solving the equation E — e2{n) — for several values of 
n. We present below analytical results for E — 62(1) = 2 for steps n > 2. Once again, we 
observe the behaviour of t and €2 around £^ = 2 for successive iterations. We set E — 2 + S, 
S being infinitesimally small, and find, by direct calculation that, for n >2 

e2{n) = ^ + fn + 0{5) (13) 

t{n) = ~+gn + 0{6) (14) 

upto leading order in 5, where /„ and Qn are respectively given by 

2^^ 2n_n 
" ~ "9~ ^ T ~ 18 



and 



i2n-2 



n 35 
3 ^ 36 



The above forms set in, as in case (i), after the first iteration, and hold perfectly well for 
n = 2, 3, 4 and 5. We now make use of the recursion relations (3) and (4) to find that the 
result is true for (n + 1)*'* stage as well. Thus, we take Eqns. (13) and (14) to represent the 
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general n behaviour of 62(71) and t{n) for any n > 2. The leading 5 - behaviour of e and i 
are obtained to be, 

1 , fn+l 



and 



1 

t = --^ + g„+i 



respectively, which lead to the equation 

^ = (15) 

i 3 3 ^ ^ 

However, it should be noted that in order that the present analysis is valid, we must have a 
finite (although large) number of iterations n such that /„ and are also finite quantities and 
are small compared to 1/5 as 5 0. The matrix elements now read (neglecting terms of the 
order of 5^ and taking the limit S 0), Pu = -2(4'* -4)/3 and, P22 = 0, Pu = -P21 = -1- 
For large (but finite) values of n one can now show, using Eqn. (5) that 



Tr,{E = 2) - 2- 



-4n 



It is quite obvious that, depending on the value of the energy, one has to select the on-site 
term and the hopping integral for the lead suitably, so that the energy does not fall beyond 
the "band" of the ordered lead or even coincide with the "band-edge" . In both cases the use 
of the expression (5) will not be meaningful. 

We have also numerically calculated T{E) for E — £2(2), E — 62(3) and E — 64(1) 
for lattices starting from TV = 3 upto N = 6. We observe a gradual attenuation in the 
value of the transmission coefficient according to our expectation. The present results thus 
provide an example of what we may call an "atypical" extended state [9] where, though the 
wavefunction displays non-zero amplitudes even at the farthest portions of an arbitrarily 
large lattice, the end-to-end transmission decays with increasing lattice size. We expect 
similar behaviour oiT{E) for other cases (with E — 622(1)) also. 
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The fixed point behaviour of T(E) : 

Before ending this section, we discuss the case where we have a fixed-point behaviour of 
the hamiltonian parameters. For this, we take a model with €22 = ^z — zt and set E = ^2 + 2t. 
All parameters then remain unaltered under RSRG and we find that for €2 = and t = 1, 
the numerical value of the transmission coefficient is given by, 

nt,) = ^-^-^ (16) 

where we have set the site energy of the lead eo, equal to zero. Naturally, we have to choose 
any suitable value for to so that the above energy remains within the "band" of the ordered 
lead. The above expression for T remains fixed for arbitrarily large versions of a diamond 
lattice. The wave function is definitely extended as t does not fiow to zero under RSRG. 

4. Conclusion 

We have presented a hierarchical lattice model where the coordination number of the 
lattice points range from 2 to 2-^"^ depending on the generation index N. In such a lattice 
there exist "extended" type of electronic states some of which have been identified and the 
corresponding eigenvalues have been calculated using renormalization group ideas. We also 
presented an exact analysis of the end-to-end transmission coefficient to reveal that the 
lattice, irrespective of its size, becomes completely transparent to an electron with energy 
E — while for other energies the transmission coefficient has a scaling behaviour. We 
obtained an exact form of the scaling for a specific energy, and the other forms can be 
obtained using the same methodology, though we did not present the other analytical results 
here. 
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Figure Captions 

Fig. 1: First tliree stages (i.e., N = 1,2,3) of the construction of a diamond hierarchical 
lattice. 

Fig. 2: Distribution of the amplitudes of an extended wavefunction at £^ = on (a) an 
N — 3 lattice and (b) a,n N — A lattice. All e^'s have been set equal to zero and t — 1. ■^j's 
take on values —1, and 1 on different sites (i). 

Fig. 3: Amplitudes of a wavefunction on (a) the basic plaquettc / which acts as a building 
block oi a,n N — A lattice and (b) an = 4 lattice for E — £2(1). All e^'s have been set 
equal to zero and t — 1. V'i's take on values —1, —1/2, 0, 1/2 and 1 on different sites (i). 
Fig. 4: Amplitude distribution on (a) an effectively N — 3 plaquette obtained by renormal- 
izing an A^ = 5 lattice twice and (b) one quarter of the original N = 5 version for E = 62(2). 
Dashed lines at the two extreme vertices indicate the presence of complementary plaqucttcs. 
All ej's have been set equal to zero and t — 1. The values of a, 6 and c are given in the text. 
Fig. 5: Amplitudes of a wavefunction for E — €4(1) on an A" = 4 lattice. All e^'s have been 
set equal to zero and t — 1. t/jiS take on values — 2-\/2, —1, 0, 1 and 2^/2 on different sites 



Fig. 6: Reduction of an n times renormalized lattice to an effective dimer. The leads are 
shown as dashed lines. 
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